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ABSTRACT 

The spectral energy balance of waves in the surf zone 
is examined with extensive measurements from the Duck94 
experiment. Cross-shore energy flux gradients are 
estimated from spectra observed with closely spaced 
pressure sensors. Nonlinear energy exchanges between 
different wave components in the spectrum are estimated 
from observed bispectra based on Boussinesq theory for 
near-resonant triad interactions. Dissipation of wave 
energy in the poorly understood breaking process is 
inferred as the residual term in the spectral energy 
balance . 

Analysis of the spectral energy balance shows that 
large decreases in energy flux observed at the dominant 
wave frequencies as waves break over a sand bar are 
closely balanced by nonlinear energy transfers to higher 
frequencies. That is, the decay of the spectral peak 
within the surf zone is a result of nonlinear energy 
transfers rather than direct dissipation. At higher 
frequencies, observed energy flux gradients are small and 
do not balance the nonlinear transfers of energy to high 
frequency components of the spectrum. This analysis 
suggests that the spectrum is saturated at high 
frequencies, and thus, the energy that cascades through 
nonlinear interactions to higher frequencies is 
dissipated in the high-frequency tail of the spectrum. 
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I. INTRODUCTION 



The lack of accurate models for the transformation of waves across the surf zone 
is the greatest deficiency is in our ability to predict surf characteristics (e.g., spilling of 
plunging breakers), near-shore circulation (e.g., longshore currents, undertow, rip 
currents), and sediment transport (e.g., beach erosion, accretion and the formation of sand 
bars). The properties of ocean surface waves on beaches are affected by changes in water 
depth, the bottom boundary layer, strong nonlinear dynamics in shallow water, and surf 
zone wave breaking. As the water depth decreases, waves refract towards the beach. As 
a result, propagation directions in the surf zone are typically close to normal incidence, 
and the shoaling evolution of waves can be reasonably well described by 1 -dimensional 
models for uni-directional waves. 

Turbulence associated with bottom friction causes a gradual decay of waves in 
shallow water, that is believed to be important to the evolution of wave spectra across 
shallow continental shelves. On beaches that typically span only 0(10) wavelengths, 
bottom damping is usually neglected. The dominant physical processes that affect wave 
transformation on beaches are believed to be nonlinear interactions and wave breaking. 
Nonlinear wave-wave interactions, weak in the open ocean, are strongly enhanced on 
beaches owing to the weakly dispersive nature of surface gravity waves in shallow water. 
Energy is continuously exchanged in near-resonant interactions between all triads of wave 
components with frequencies fi,f 2 ,fj, that obey the interaction rule fj ± fj ± fs = 0 {Freilich 
and Guza, 1984). These so-called triad interactions not only cause a broadening of the 
wave spectrum (i.e., energy transfers from the dominant waves to higher-and lower 
frequencies), but are also responsible for the dramatic transformation of wave shapes from 
nearly symmetric sinusoidal deep water swells to the pitched-forward wave shape of near 
breaking waves on beaches {Elgar and Guza, 1985a). 

Models for wave shoaling transformation outside the surf zone are usually based 
on nondissipative Boussinesq equations for weakly nonlinear, weakly dispersive waves 
{Peregrine, 1961). These models include time-domain formulations {Madsen et al., 1991), 
discrete frequency domain formulations {Freilich and Guza, 1984; Liu et al., 1985) and 
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stochastic formulations that predict the evolution of continuous wave spectra {Berbers and 
Burton, 1997; Norheim et a/., 1998). Extensive comparisons to field data confirm that 
Boussinesq models predict accurately the shoaling evolutions of waves outside the surf 
zone {Freilich and Guza, 1984; Elgar and Guza 1985b; Norheim et al., 1998). 

Early surf zone models [Collins, 1970; Battjes, 1972; Goda 1975] simply 
parameterized the heights of breaking waves in the terms of the local water depth. Battjes 
and Janssen [1978] introduced a frequency- integrated energy balance model with the 
dissipation sink term based on the energy loss in a turbulent bore and a simple breaking 
criterium. Thornton and Guza [1983] refined this model by incorporating realistic wave 
height statistics. Extensive comparisons with laboratory and field data for a range of 
beach profiles and incident wave conditions shows that the energy balance models of 
Battjes and Janssen and Thornton and Guza predict accurately the variations of wave 
heights in the surf zone {Battjes and Janssen, 1978; Thornton and Guza, 1983; Battjes 
and Stive, 1985). Although wave height decay in the surf zone is predicted accurately by 
these bulk models for total wave energy, the associated dramatic changes in the wave 
spectrum resulting from nonlinear interactions are not described by these models, and 
therefore the characteristics of breaking waves (e.g., spilling vs plunging breakers; narrow 
vs broad spectra) are not predicted. 

Recently, Boussinesq models, that rigorously account for nonlinear interactions, 
were extended into the surf zone by including a heuristic dissipation term in the shoaling 
evolution equations {Mase and Kirby, 1992; Kaihatu and Kirby, 1995; Eldeberky and 
Battjes 1995; Chen et al., 1997). As virtually no information exists on the spectral 
characteristics of dissipation rates in the surf zone, simple parameterizations were used. 
Mase and Kirby, (1992) and Kaihatu and Kirby, (1995) used a frequency dependent 
dissipation function that combines a constant term that reduces spectral levels uniformly 
in proportion to the energy at each frequency with a frequency- weighted (f^) term that 
dissipates the higher frequency components, using an adjustable coefficient to control the 
relative importance of the two terms. Eldeberky and Battjes (1995) used a frequency 
independent formulation, in which dissipation is proportional to the energy at each 
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frequency. Chen et ai, [1997] showed that the evolution of wave shapes across the surf 
zone is predicted more accurately by models with a dissipation term that increases with 
frequency than models using a frequency independent dissipation term. 

To obtain a better imderstanding of the spectral characteristics of energy 
dissipation within the surf zone, the spectral energy balance of shoaling and breaking 
waves is examined in this study using extensive field observations from a natural sandy 
ocean beach. Neglecting alongshore depth variations, directional spreading of waves and 
wave reflection from shore, the energy balance can be expressed as 



where the left-hand side is the cross-shore gradient of the energy flux F(f). The source 
term S„^ represents the net nonlinear transfer of energy to waves with frequency / owing 
to triad interactions. The other source term Sj^(f) represents energy dissipation during 
wave breaking. The energy flux gradients are evaluated from closely spaced instruments. 

The energy transfers in near-resonant triad interactions are a function of the 
phase-relationships between the three interacting components, and thus depend on higher- 
order statistics of the wave field. A stochastic formulation of the Boussinesq equations 
yields the following expression for S„^ (Herbers and Burton,l997; Norheim et al., 1998) 




( 1 ) 






(2a) 




(2b) 



0 



0 



3 



where IM { } indicates the imaginary part. The third-order bispectrum B(f,,f 2 ) describes 
the degree of coupling and the phase relationship in triads of wave components with 
frequency f,,f 2 and fi+f 2 (Hasselmann et. ai, 1963). In shallow water, the bispectrum 
evolves strongly and describes in a statistical sense the shapes of shoaling waves (e.g., 
Elgar and Guza, 1985b; Elgar et ai, 1990a). The nonlinear source term 5„,(/) in the 
energy balance involves an integral b(f) of over aU triads that include a component 

with frequency f. 

In the present study the terms in the energy balance equation (1) were examined 
with measurements from a cross-shore array of pressure sensors deployed near Duck N.C., 
in 1994. The energy flux gradient F,(/) and nonlinear source term 5„/(/) were estimated 
from measured spectra and bispectra at 14 cross-shore positions. The associated 
dissipation rate SJIJ) was then inferred as the residual term in equation 1. The field data 
and analysis are described in section 2. 

Example observations of the energy balance in a range of conditions are presented 
in section 3, followed by a statistical analysis of the entire 2 months observations in 
section 4. In the surf zone, large negative values of FJf) in the energetic part of the 
spectrum are shown to be approximately balanced by 5„,(/). These results show that the 
energy of the dominant waves is not directly dissipated in the surf zone, but transferred 
to higher frequencies through nonlinear interactions where the energy is presumably 
dissipated. The results are summarized in section 5. 
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II. FIELD DATA AND ANALYSIS 



Detailed field measurements of the evolution of shoaling and breaking waves on 
a natural ocean beach were acquired during the Duck 94 experiment at the U.S. Army 
Corps of Engineer Field Research Facility located near Duck, N.C., on a relatively straight 
barrier island exposed to the open Atlantic ocean. A dense cross shore transect of 14 
pressure sensors was deployed, extending from the shore line to about 5 m depth (350 m 
from shore)(Figure 1). Collocated with each pressure sensor was a sonar altimeter that 
provided accurate estimates of the location of the sea floor. The sample frequency of all 
sensors was 2 Hz. High quality data were collected nearly continuously during the 
months of September and October 1994 (Elgar et al, 1997). 

The closely spaced pressure sensors allowed for direct estimates of the advection 
term F^(f) and the nonlinear source term in the energy balance equation (1). A 
surface elevation spectrum E(f) and bispectrum B(f,,f 2 ) with resolution of .0078 Hz were 
estimated from each hour-long pressure record using the shallow water approximation. 

The spectral energy balance was evaluated at the midway point between adjacent 
instruments using a simple finite difference approximation. The energy flux gradient F^(f) 
was approximated by (using the shallow water approximation of the group speed): 



1 
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X -X 
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(3) 



where Xjpc 2 (>Xj) are the cross-shore positions of a pair of adjacent sensors, hj and hj are 
the corresponding water depths, E,(f), and E 2 (f) are the corresponding spectra. 

The nonlinear energy transfer term (Sqi) at the same midway point was 
approximated by: 






( 4 ) 
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with b,(f) and b 2 (f) the bispectral integral b(f)(Eq.2b) evaluated at the instrument locations 
X, and Xj. 

Significant weather changes throughout the Duck 94 experiment provided ample 
sampling of a multitude of surf conditions including two major Noreaster storms with 
maximum significant wave heights of 2.5 m and 3.8 m (Figure 2) and surf zones 
extending across the entire instrumented transect. The data set also includes calm periods 
with significant wave heights as small as .2 m and virtually no breaking across the 
transect. Observed wave spectra include narrow spectra of long period, remotely 
generated swells, broad spectra of locally generated seas and mixed swell-sea spectra. 
The bottom profile during the experiment consisted of a mild offshore slope with a well 
developed bar approximately 100 m from shore followed by a trough and steeper 
foreshore. (Figure 1). Fluctuations of bottom topography in the surf zone were minor 
throughout the data collection period with the exception of the passage of the largest 
storm that moved the bar 80 meters seaward. 
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III. OBSERVED SPECTRAL ENERGY BALANCE 



The spectral energy balance of shoaling and breaking waves observed for a wide 
range of conditions is illustrated here with four case study examples (Figure 2), each 
based on a 1-hour long data record. For each case, the evolution of the spectrum E(f) 
across the beach is shown in Figure 3, and a representative video image of the sea surface 
in the instrumental area is shown in Figure 4. In case I, small amplitude swell propagated 
across the instnomented transect with little breaking (Figure 4). The evolution of the 
narrow spectrum shows the expected growth of harmonic peaks. Breaking and associated 
energy losses occurred primarily on the beach face shoreward of instnunent array (Figure 
4). 

Case studies 11 and III are from the two major storm that passed through the 
region during the experiment. In both cases the irregular surf zone with spilling breakers 
extended across the entire 350m-long transect (Figure 4). The evolution of the spectra of 
these broad wind seas across the surf zone is similar with a strong, nearly uniform decay 
of spectral levels in the energetic part of the spectrum (Figure 3). In case IV, during the 
waning of the second Noreaster, moderately energetic swell was observed. This case is 
characterized by a narrow surf zone confined to the iimer part of the instrument array, 
with some intermittent breaking occurring farther offshore on the sand bar (Figure 4). The 
spectra observed in case IV indicate a pronounced second harmonic peak outside the surf 
zone. Inside the surf zone both the primary spectral peak and the second harmonic peak 
are attenuated. 

To examine the role of nonlinear interactions and dissipation processes in the 
observed spectral evolution of these four case studies, estimates of the energy flux 
gradient F,(f) (Eq. 3) and the nonlinear transfer S„j(f) (Eq. 4) are compared in Figmes 5-8. 
The dissipation term Sds(f)in the energy balance can be inferred as the difference between 
the estimated F;j(f) and S^(f) (Eq. 1). In the non-breaking wave Case I (Figure 5) the 
observed values of S^j and F^ are approximately equal at all locations except for the 
innermost two, indicating that the observed spectral changes (Figure 3) are the result of 
nonlinear energy exchanges, consistent with Boussinesq model results for similar 



7 



conditions presented in Elgar et al, (1997) and Norheim et al, (1998). The positive 
values of S„| and in the range 0. 1-0.4 Hz indicate a transfer of energy from the 
dominant 0.07 Hz swell to high-frequency harmonic components, consistent with the 
Boussinesq model predictions of Elgar et al. (1997) and Norheim et al, (1998). The 
largest transfers occur slightly inshore of the bar crest (H,I) and on the beach face (L). 
Comparisons of F, and at low frequencies (<0.1 Hz) are not shown because partial 
reflection of the small amplitude swell from the beach face (Elgar et al., 1997), caused 
large errors in the energy flux estimates (i.e., standing wave nodes and antinodes). 
Reflection effects are apparent in the large fluctuations of F, observed at the innermost 
stations K and L (Figure 5). 

In case II with energetic seas breaking across the entire transect (Figure 6), the S^j 
estimates are characterized by a negative lobe centered at the spectral peak frequency, and 
a broader positive lobe at higher frequencies. The integral of S,j(f) over all frequencies 
is zero because the triad interactions conserve energy. The Snj estimates show that energy 
is primarily transferred from the energetic part of the spectrum to higher frequencies. The 
observed energy transfer increases from the outer most station C (A and B were not 
operational during this event) to a maximum value on the sand bar crest (F) decreasing 
to small values near the shoreline (i.e., where most of the energy has been dissipated). 
The estimated energy flux gradients F^ show a negative lobe in the energetic part of the 
spectrum that nearly equals the negative lobe in the estimates. The surprisingly close 
balance of F, and demonstrates that the large energy losses observed in the energetic 
part of the spectrum (Figure 3) are accounted for by nonlinear transfers to higher 
frequencies, and thus no significant dissipation takes place in the energetic part of the 
spectrum. 

The observed energy flux gradients and nonlinear transfers do not balance at 
higher frequencies. Whereas Sni(f) estimates show a broad positive lobe above the spectral 
peak frequency, the associated Fj(f) estimates are close to zero. These results suggest that 
the energy transferred through nonlinear interactions to higher frequencies cannot be 
absorbed in the high-frequency tail of the spectrum and is presumably dissipated. 
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qualitatively consistent with the saturation of wave spectra in the surf zone observed by 
Thornton (1977, 1979). In the most energetic Case HI, similar fuUy developed seas with 
a lower peak frequency (.09 Hz) were observed (Figures 3,4). In this severe storm, several 
instruments failed and thus the energy balance could not be evaluated on the seaward side 
of the bar crest. Nevertheless the estimates of Sni(f)and F,(f) (Figure 7) are remarkably 
similar to the case II estimates, with a close balance between energy losses in the 
energetic part of the spectrum (the negative lobe of F^(f))and energy transfers to higher 
frequencies (the negative lobe of Snj(f)), and an apparent saturation of the spectrum at 
high frequencies (S^j > 0 and F, ~ 0). As in Case II the largest energy transfers are 
observed at station F on the bar crest. At the innermost stations K and L, fluctuations in 
F^(f) at lower frequencies indicate that wave reflection from the steep beach face (i.e., 
nodes and antinodes) caused significant errors in the energy flux estimates. 

A few days later (Case IV) the storm has moved offshore and swell with a peak 
frequency of 0.09 Hz (Figure 3) propagated though the array. Seaward of the bar crest 
at stations A and B energy is transferred from the spectral peak to it’s secondary harmonic 
and higher frequencies (Figure 8). The negative values of F, observed at station B do not 
balance the positive Snj estimates, suggesting some energy is dissipated. No data was 
collected on the bar crest but the large reduction in spectral levels between stations B and 
F (Figure 3) indicates strong dissipation takes place in intermittent wave breaking on the 
bar crest. Between the bar and the beach face at stations F-J, nonlinear energy transfers 
are approximately balanced by the estimated energy flux gradients. Whereas at station F, 
energy is transferred from the primary peak to the 2“** harmonic, further inshore at stations 
I and J energy is transferred from the second harmonic to higher frequencies. Finally, at 
the shallowest stations K and L located well inside the surf zone, the estimates show 
large transfers of energy from the energetic part of the spectnmi to higher frequencies. 
As in cases II and III, the F, estimates closely balance in the energetic part of the 
spectrum, and F, « 0 at high frequencies, indicating that the decay of the energetic part 
of the spectrum results from nonlinear transfers to the saturated high-frequency tail of the 
spectrum where the energy is dissipated. 
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IV. DISCUSSION 



The nonlinear transfer of energy to or from wave components with frequency f is 
controlled by the phase 0(f) of the bispectrum integral b(f) (Eq. 2). In intermediate water 
depths (kh=0(l)) the theoretical value of this so called biphase is 0“ or 180® (Hasselmarm 
et al., 1963). For 0 = 0 or 180®, the imaginary part of b(f) = 0 and thus no transfers of 
energy take place. At these depths the bispectrum describes the nonresonant coupling 
between the primary wind waves and bound secondary waves. In shallow water (kh«l) 
biphases evolve allowing for energy transfers in near-resonant triad interactions. Negative 
biphase values correspond to a loss of energy for components with frequency f, and 
positive values (O<0(f)<18O®) indicate energy gains. The biphase 0 also characterizes the 
shape of waves that typically evolve from symmetric profiles in deep water (0=0 or 180®) 
to pitched-forward crests prior to breaking in shallow water (0= ± 90®) (Elgar and Guza, 
1985). 

Biphase estimates are shown in Figure 9 for all four case studies at locations: 
seaward of the bar, near the bar crest, and on the beach face. Seaward of the bar, the 
biphases are close to the intermediate depth values 0® or 1 80®. Case I shows more scatter 
probably because the nonlinearity is weak in this case and thus the bispectral estimates 
have greater statistical uncertainty. Near the bar crest, the bisphases have evolved in 
cases I-in to about 135® near the peak frequency and to 30-90® at higher frequencies. The 
strongest asymmetry is observed in the non-breaking Case I where 0 is close to 90® at the 
second harmonic peak frequency. Interestingly in cases I-III biphases observed on the 
beach face show a relaxation towards the values 0 and 180® for symmetric wave profiles. 
On the other hand, in case IV biphases observed inshore of the bar crest (where waves 
reform after breaking on the bar) are close to 0 and 180®, and strong asymmetry does not 
develop imtil waves break again on the beach face. These results indicate that 
asymmetries in wave profiles are most pronounced during or before the onset of wave 
breaking and disappear in the inner region of wide surf zones. 

The energy balance observed in case studies II-IV with energetic breaking waves 
(i.e., the surf zone extends seaward of the sand bar) is qualitatively similar with nonlinear 
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interactions transferring energy from the energetic part of the spectrum to higher 
frequencies where the energy is dissipated. This surf zone balance is shown in Figures 10 
and 1 1 to hold for all observations with an energetic wide surf zone during the two month 
long experiment. Only 77 low-tide observations with incident wave variances measured 
in 8 m depth > 2000 cm^ (i.e., significant wave heights > 1.8 m) are included to restrict 
the comparisons to wide surf zones that extend across the entire instrumented transect. 
Estimates of F^(fp) and Snj(fp) at the spectral peak frequency (fp) are compared at locations 
A-L(midway between adjacent instruments) in Figure 10. 

Between the seaward end of the transect and slightly inshore of the bar crest 
(locations A-1), S„](fp) and F,(fp) are both negative and approximately equal. Correlation 
coefficients between Snj and F, vary between .76-. 96 with least-squares fit lines that are 
close to (within 5-30%) a one-to-one correspondence. These estimates demonstrate that 
the large energy losses observed in the surf zone at the spectral peak frequency are 
accounted for by nonlinear transfers to higher frequencies rather than direct dissipation 
of wave energy in the energetic part of the spectrum. 

Qose to shore (locations J-L) the Snj(fp) estimates are small and in poor agreement 
with Fj(fp) estimates. At these inner surf zone locations, energy levels may be reduced 
enough that significant reflected components from the beach (e.g., an antinode near the 
shoreline) contribute large errors to the energy flux gradient F,j(fp) estimates (e.g., the 
large positive F^(fp) estimates at L). 

To characterize the energy balance at higher frequencies, the S„,(f) and F,^(f) 
integrated over the frequency range of .3-.5 Hz are compared in Figure 11. Whereas the 
Sn] estimates show consistently positive values (i.e., energy transfers to high frequencies, 
F,j estimates are comparatively small. These results confirm that the high-frequency tail 
of the spectrum is saturated in the surf zone and the energy transferred to high 
frequencies in nonlinear interactions is irmnediately dissipated. The largest energy 
transfers (and hence dissipation rates) are observed at stations F-L, near the crest of the 
sand bar, and L on the beach face where the strongest breaking activity is expected. 
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V. SUMMARY AND CONCLUSIONS 



The spectral energy balance of shoaling and breaking waves was examined with 
two months of field observations from a 350m-long cross-shore transect of 14 pressure 
sensors deployed on a barred ocean beach near Duck, NC as part of the Duck 94 
Experiment {Elgar et al., 1997). The energy flux gradient term Fj,(f) and the nonlinear 
source term S„,(f) in the energy balance equation were estimated at 13 cross shore 
positions from approximately 1400 hour-long pressure records. The energy flux gradient 
Fjj(f) was estimated from the measured changes in the wave spectrum E(f) between two 
adjacent instruments. The net nonlinear energy transfer Snj(f) resulting from near-resonant 
triad interactions was estimated from the observed wave bispectrum B(f,,f 2 ) based on a 
stochastic Boussinesq model {Herbers and Burton, 1997). 

In low-energy conditions, when waves propagate across the instrumented transect 
with little or no breaking, the observed small nonlinear energy transfers to higher 
frequencies are approximately balanced by the growth of the spectrum at high frequencies, 
consistent with earlier studies Elgar et al.{\991), Norheim et a/.(1998). In high-energy 
conditions with waves breaking across the entire transect, the Sn,(f) estimates consistently 
show a large negative lobe at the spectral peak frequency and a broad positive lobe at 
higher frequencies, consistent with a strong nonlinear transfer of energy from the spectral 
peak to a broad range of higher frequencies. The largest energy exchanges are observed 
on the crest of the sand bar, where wave breaking is most intense. The observed negative 
energy flux gradients in the energetic part of the spectrum closely balance the nonlinear 
transfers. Theses observations show that the large energy losses in the spectral peak of 
breaking waves are accounted for by nonlinear transfers to higher frequencies, and thus 
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no significant dissipation occurs in the energetic part of the spectrum. At higher 
frequencies the observed are small and do not balance the positive Sn, values. These 
results support earlier observations {Thornton, 1977, 1979) that wave spectra in the surf 
zone are saturated at high frequencies (i.e., energy transferred through nonlinear 
interactions cannot be absorbed in the high-frequency tail of the spectrum and is 
presumably dissipated). 

In conclusion, the present observations of the spectral energy balance in the surf 
zone show that near-resonant triad interactions transfer energy from the energetic part of 
the wave spectnun to higher frequencies where the energy is dissipated. 
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APPENDIX 



Figure 1. Instrument locations in the Duck94 experiment. Squares indicate the 
14 cross-shore positions of a co-located pressure sensor and sonar altimeter. A 
(outermost) through M (innermost) indicate the midway points between adjacent 
instruments. 

Figure 2. Variability of incident wave variances estimated in the frequency range 
0.05-0.25 Hz from bottom pressure measiuements in 8 meter depth (instnunents deployed 
and maintained by the U.S. Army Corps of Engineers) using a linear theory depth 
correction. Hourly variance estimates are shown versus time for the two months of 
continuous data collected. Dotted vertical lines denote the four case studies analyzed in 
section 3. 

Figure 3. Observed shoaling evolution of wave spectra in the four case studies. 
In each case spectra at foiu locations spaiming the transect (Figure 1) are shown. 

Figure 4. Representative surf zone photographs (courtesy of T.C. Lippmarm) of 
four case studies obtained with video cameras mounted on a nearby tower about 50 m 
above the sea surface. The top, middle-left and lower-left photographs all cover the same 
area including the inner half of the instrumented transect (solid line). The middle-right 
and bottom-right photos provide simultaneous, slightly overlapping (!=20m) images that 
cover the outer half of the instrumented transect (sohd line). The outward looking camera 
was not operational during cases I and II. 

Figure 5. Estimates of the energy flux gradient Fj^(f)(dot-dash) and the nonlinear 
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transferred energy Stj(f) (solid) versus frequency f at 10 cross-shore locations in case 
study I. The beach profile and the cross-shore locations where the energy balance is 
evaluated are indicated in the bottom right panel. Arrows indicate the spectral peak 
frequency. 

Figure 6. Estimates of Fj(f) and Sn](f) in Case study II at 10 cross -shore locations 
(same format as figure 5). 

Figure 7. Estimates of F,(f) and Sn,(f) in Case study in at 8 cross-shore locations 
(same format as figure 5). 

Figure 8. Estimates of F^(f) and Sni(f) in Case study IV at 9 cross-shore locations 
(same format as figure 5). 

Figure 9. Biphase 0 versus frequency f (Hz) seaward of the bar, near the bar 
crest, and on the beach face. From top to bottom; case I at locations C,G,L, case II at 
C,G,L, case III at A,F,L, and case IV at A,F,L,. 

Figure 10. Nonlinear energy transfer versus the energy flux gradient F, at the 
spectral peak frequency in energetic surf zone conditions. Each scatter diagram shows the 
comparison at one cross-shore location (comparisons at locations D and E where only few 
observations were collected are not shown). Each dot represents an estimate based on a 
1-hour-long data record. A linear regression line (dash-dot) and correlation coefficient 
between F^ and are indicated in each scatter diagram. For reference, a representative 
beach profile with the instrument locations is shown in the lower right panel. 

Figure 11. Nonlinear energy transfer S^i versus the energy flux gradient F,^, both 
integrated over the high-frequency range 0.3-0.5 Hz, in energetic surf zone conditions 
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(same format as figure 10). 
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